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ABSTRACT 

About one year ago, it was speculated that decaying or annihilating Light Dark Mat- 
ter (LDM) particles could explain the flux and extension of the 511 keV line emission 
in the galactic centre. Here we present a thorough comparison between theoretical 
expectations of the galactic positron distribution within the LDM scenario and obser- 
vational data from INTEGRAL/SPI. Unlike previous analyses, there is now enough 
statistical evidence to put tight constraints on the shape of the dark matter halo 
of our galaxy, if the galactic positrons originate from dark matter. For annihilating 
candidates, the best fit to the observed 511 keV emission is provided by a radial 
density profile with inner logarithmic slope 7 = 1.03 ± 0.04. In contrast, decaying 
dark matter requires a much steeper density profile, 7 > 1.5, rather disfavoured by 
both observations and numerical simulations. Within the annihilating LDM scenario, 
a velocity-independent cross-section would be consistent with the observational data 
while a cross-section purely proportional to can be rejected at a high confidence 
level. Assuming the most simplistic model where the galactic positrons are produced as 
primaries, we show that the LDM candidate should be a scalar rather than a spin-1/2 
particle and obtain a very stringent constraint on the value of the positron production 
cross-section to explain the 511 keV emission. One consequence is that the value of the 
fine structure constant a should differ from that recommended in the CODATA. This 
is a very strong test for the LDM scenario and an additional motivation in favour of 
experiments measuring a directly. Our results finally indicate that an accurate mea- 
surement of the shape of the dark halo profile could have a tremendous impact on the 
determination of the origin of the 511 keV line and vice versa. 
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1 INTRODUCTION 

An emission line at 511 k eV was detected at the galactic 
centre three decades ago i Johnson et alJ^1972^ . Its identi- 
fication as an electron-positron annihilation line followed 
as soon as high-resolu tion spectrometers became available 
iLeventhal et alJl978^ , but the origin of low-energy galactic 
positrons is still a matter of heated debate. The latest ob- 
servations of the annihilation emission have been performed 
by the SPI spectrometer aboard the INTEGRAL^ satellite. 
A total flux of ~ 10""^ photons s~^ cm~^ was measured, 

* E-mail: yago@head.cfa.harvard.edu 

^ INTEGRAL (International Gamma Ray Laboratory) is an 
ESA's gamma ray observatory launched in October 2002. 



in agreement with previous estimates. The morphology of 
the galactic bulge emission could be fit by a Gaussian with 
~10 °FWHM. A disc component was recently detected by 
INTEGRAL/SPI (Knodlseder et al. 2005) but this emission 
can be attributed to the (i^ decay of the radioactive species 
^®A1 and ''''Ti, which are produced by massive stars in the 
disc. 

Several astrophysical sources have been proposed in the 
literature to explain the low-energy positrons from the bulge, 
such as radioactive nuclei expelled by stars (supernovae, hy- 
pernovae, novae, Wolf-Rayet stars and red giants) and col- 
lapsed objects (neutron s tars or black holes). Ne vertheless, 
most of these sources (see lKnodlseder et all2005l . and refer- 
ences therein) cannot account for the observed morphology. 
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due to the large bulge-to-disc ratio of the emission, which 
suggests an old stellar population origin, unless rather elab- 
orate mechanisms (e.g. jets, propagation) are invoked. 

On the other hand, the presence of low-energy positrons 
could be exp l ained by Dark Mat ter (DM) annihi l ations 
jBoehm et alJ 1200411 or decays jHooper fc Wand 12004 
iPicciotto fc Pospelovlliooil) . The present paper focuses on 
such scenarios, which require light dark matter particles (i.e. 
with a mass mdm J; 100 MeV, depending on their exact na- 
ture) in order to reproduce the observational data. 

The smallness of the DM mass might appear surpris- 
ing to many. Indeed, not so long ago, most of the com- 
munity thought that annihilating DM particles should be 
heavier than a few GeV because of the Lee- Weinberg limit 
jLee fc Weinberslll977ll . which states that if DM is a sta- 
ble fermion coupled to heavy particles (such as the Z and 
W gauge bosons) then its mass should exceed that of the 
proton; otherwise it would overdose the universe. There are 
possible ways to evade the Lee- Weinberg limit though but 
the fact that theoretically motivated DM candidates, such 
as the lightest neutralino, were naturally very heavy did not 
encourage the community to investigate the lighter range. 

The window for Light Dark Matter (LDM) particles 
suddenly opened when it was realized that scalar candidates 
with a mass from a few MeV to a few GeV, coupled to heavy 
fermions [F) or to light neutral particles (neutral gauge 
bosons Z' , somewhat analogous to the Z gauge boson), could 
also yield the observed relic density. However, the introduc- 
tion of LDM particles immediately faces an embarassing 
problem: their annihilations are expected to produce too 
many low-energy gamma rays in our galaxy, compared to 
what has been observed. To be on the safe side, the present- 
day annihilation cross-section must be about five orders of 
magnitude (times m|j(,v, where rriMoV = jridmC'^/l McV) 
small er than it was in the primordial universe (|Bochm ct al. 
120041) . 

Such a condition can be easily satisfied if the channel 
associated to the exchange of heavy fermions is suppressed 
with respect to the one due to the new gauge boson. Indeed, 
the annihilation cross-section associated with a Z' exchange 
is proportional to the square of the DM velocity, which - in 
the Milky Way - is at least two or three orders of magni- 
tude smaller than in the primordial universe (i.e. before DM 
became non-relativistic) . Hence a velocity-dependent cross- 
section can satisfy both the relic density criterion and the 
gamma ray constraint. In contrast, the cross-section arising 
from heavy-fermion exchange does not depend at all on ve- 
locity; it remains constant at any epoch. The LDM scenario 
is thus viable if - as initially proposed - the F exchange 
is at least five orders of magnitude (times m|ij,v) smaller 
than the Z' cross-section at early times. Note, though, that 
the contribution of the fermion exchange to the total anni- 
hilation cross-section could become dominant as velocities 
become non-relativistic. 

Although the parameter space allowed for LDM is still 
quite broad, it is encouraging to find out that this model, 
originally built to be invisible (more precisely, to escape 
the low-energy gamma ray constraint) is able to explain 
the observed properties of the 511 keV line without adding 
any new components nor changing the mass range initially 
proposed. Moreover, the existence of LDM particles might 
also explain several particle physics measurements (see e.g. 




Figure 1. Dark matter annihilation through the exchange of a 
charged heavy fermion F (left) or a neutral light gauge boson Z' 
(right). 

iBoehm fc Ascasiba dliooi) . The detailed analysis of the 511 
keV line presented here provides new - and relatively tight 
- constraints on some of the LDM parameters, greatly en- 
hancing the predictive capability of the model in other fields. 

Conversely, independent confirmation of the LDM sce- 
nario from particle physics experiments would have impor- 
tant astrophysical consequences. In fact, if dark matter anni- 
hilation (light or otherwise) turns out to be the main source 
of galactic low-energy positrons, the observed morphology 
of the 511 keV emission line would constitute an excellent 
probe of the shape of the Milky Way dark matter halo, what- 
ever the exact nature of dark matter particles might be. The 
present study shows that several robust conclusions can al- 
ready be derived from curre nt INTEGRAL/S PI data. 

Previous analyses (e.g. iJean et'al] 120031) have shown 
that a point source can be ruled out a high confidence level. 
If positrons cannot travel a long distance before annihilat- 
ing, that would seem to suggest that the galactic halo cannot 
be too 'cuspy'. On the other hand, a very flat profile would 
not match the observed morphology of the emission either, 
because its FWHM tells us that most of the positrons are 
ge nerated within 1 kpc from the galactic centre. According 
to lBoehm fc Ascasibar! i2004l) . the observed morphology is 
well described if the dark matter halo of the Milky Way fol- 
lows a Navarro ct al. ( 1997) profile. In this paper, we will 
attempt to constrain both the nature of dark matter and its 
distribution within the Milky Way halo as independently as 
possible. 

Section 121 focuses on the main features of our dark mat- 
ter model. The description of the Milky Way halo is dis- 
cussed in Section |3 Section 21 is devoted to the compari- 
son between the theoretical positron distribution and the 
observed fiux. The results of our likelihood analysis are re- 
ported in Section 1^1 while Section |S| discusses additional is- 
sues, indirectly related to the present work. Finally, our main 
conclusions are briefiy summarized in Section Q 



2 DARK MATTER CHARACTERISTICS 

Our study is based on the assumption that most galactic 
positrons originate from the decays or annihilations of LDM 
particles. The number density of positrons produced per unit 
of time is then dictated by the number density of dark mat- 
ter particles, Wdm, times their decay/annihilation rate into 
a pair e^e~ . The latter are given by Fd or Fa = {(JVr)ndui* 
respectively, where {avr) is the thermal average of the an- 
nihilation cross-section times the DM relative velocity. The 
two annihilation channels mentioned in the introduction are 
represented in Figure Q 

Since the energy of the particles in the final state is im- 
posed by kinematics, the only quantities that may appear in 
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the formula of the annihilation cross-section are the masses 
of all particles involved, their couplings and the energy -Edm 
of the initial state. Since DM particles are nowadays non- 
relativistic, £dm ~ mdm c? +m-dm «^/2, so the average cross- 
section can be re-written as {avr) ~ a -I- hv^ -\- 0{v'^), where 
both V and Vr are expressed in units of the speed of light, c. 

In our case, the _F-exchange cross-section has both an a- 
and a fe-term; both of the same order of magnitude. In con- 
trast, the Z' exchange gives rise to a pure velocity-dependent 
cross-section (i.e. a = Q). A combination of both ingredi- 
ents [Z' and heavy fermions) thus provides a total annihi- 
lation cross-section with a ^ b ^ Q. Moreover, the cross- 
section through the Z' exchange depends on rridm, while 
the cross-section associated with the exchange of heavy par- 
ticles is proportional to l/m% or m^,-„/mf-, depending on 
whether DM is a spin-0 or spin-1/2 particle, respectively 
(Bochm & Favct 2004). 

There are two independent constraints on the values 
of a and b (or more precisely on {aVr}). On one hand, the 
relic density criterion imposes that the annihilation cross- 
section at the time of the chemical decoupling (i.e. when 
T ~ mdm and Udm ~ c) is about (avr) ~ 10~^^ cm^ 
so that fldmh^ ~ 0.1 nowadays. On the other hand, the 
amount of low-energy gamma rays produced at the present 
day in the centre of the Milky Way should n ot exceed the 
value s observed by COMPTEL and EGRET iStrong et alJ 
I2OOOI) . An estimate of the gamma ray flux generated by 
LDM annihilations has been computed by iBoehm et alJ 
l|2004) in terms of {ovr), yielding^ an upper limit {uVr) < 
10~^^mMev cm^ s~^. Combining both constraints, we ob- 
tain 

(cr«Oprim « a -f &/9 ~ 10"^'' cm^s"^ (1) 

in the primordial universe (with iidm ~ c/3) and 

{cr«r)now ~ a -f -Uq b < lO^^^mMev cm^s""^ (2) 

in our Galaxy, where vo ~ 10~'^c and assuming that ev- 
ery DM annihilation into electron-positron instantaneously 
produces two photons with the maximal energy. Given that 
the photon production is expected to be through the final 
state radiation (i.e. a radiative correction to the annihilation 
process), we are clearly overestimating the gamma ray flux, 
making the constraint more stringent than it should be (see 
the discussion in Section 1^!^ . 

In any case, it seems clear that, for mdm ^ 100 MeV, 
the annihilation cross-section at the present time must be 
suppressed with respect to its primordial value in order to 
satisfy both the relic density Q and the gamma ray ^ 
conditions, which lead lead to a < 10"''^ m^cv ^ni^ s~^ and 
b ~ 10~^^ cm^ s~^. A velocity-dependent cross-section is 
thus necessary for LDM candidates below ~ 100 MeV. 

The gamma-ray constraint depends nevertheless on the 
shape of the dark halo profile. In particular, equation ^ 
above was derived from a NEW profile. Taking a fiat- 
ter halo lessens these constraints and one finds that a 



The flux obtained by iBoehm et alj l200'd) was actually overes- 
timated by at least a factor four, since it was assumed that the 
number densities of dark matter particles and anti-particles were 
given by ridm = "dm* = Pdm/^dm instead of ridm = "dm* = 

Pdm/(2m.dm)- 



velocity-independent cross-section can match both the relic 
density and the gamm a ray criteria if mdm > 20 MeV 
jAhn fc KomatsulboO^ . Hence, if the dark halo radial den- 
siy profile of the Milky Way turns out to be flat, it would 
be possible to set & = and get rid of the Z' . 

In the following, the constants a and b will be normal- 
ized to 10"^" cm^ s~^, yielding the notation a26 and 626, 
respectively. For decaying LDM part icles, the de cay rate Fd 
will also be normalized to 10^ -'^^ s"^ iHooper fc Wane 200l : 
iPicciotto fc Pospelovllioosl) yielding the notation r26. 



3 THE MILKY WAY DARK HALO 

Pioneering analytical studies based on the sp herical collapse 
formahsm jGunn fc Gotlill972l: iGunnlllQTTt) predicted that 
dark matter haloes ought to be described by a single power- 
law density profile, p(r) oc r \ w i th 7 ranging from 2 to 2.25 
(e.g. iFiUmore fc Goldreichl Il984l: lBertschingeilll985l) . Such 
result appeared quite compelling at the time, as it could 
straightforwardly explain the flatness of the rotation curves 
observed in spiral galaxies. 

Later work based on numerical N-body simulations 
showed that the density profile was shallower than isother- 
mal (7 < 2) as r — » 0, and steeper (7 « 3) as r ^ 00. Typ- 
ical values measured for the asymptotic logarithmic slope 
at the centre range fro m 7=1 fcj avarro ct al. 1997, here- 
after NEW) to 7 = 1.5 jMoore et alJll999l . hereaft^er M99). 
The very existence of an asymptotic beh aviour has recently 
been questioned by several studies (e.g . IPower et alj l2003l : 
iHavashi et"al]l2004l: iNavarro et alj|2004t) . in which the den- 
sity profiles found in high-resolution simulations are re- 
ported to become progressively shallower inwards. 

Despite the significant uncertainty on the shape of the 
density profile near the centre^, there is general agreement 
in that the dark matter distribution within a spherically 
symmetric halo can be well fitted by a 'universal' function 
with a small number of free parameters, and that the same 
functional form is valid for a broad range of halo masses 
and underlying cosmologies. Most of the analytical formulae 
proposed in the literature can be cast in the form 



(r/ro)^[l + (r/ro)"l'''-^'/"' 

where po and ro are a characteristic density and radius of 
the halo, 7 is the asymptotic logarithmic slope at the centre, 
P is the slope as r ^ 00 and a controls the exact shape of 
the profile in the intermediate regions around ro. 

Many different sets of values have been suggested for 
these parameters. The most notable discrepancies concern 
the value of 7. In particular, observed rotation curves of 
dwarf spiral and low surface brightness galaxies tend to 
favour fiat profiles (7^0), which has bee n often signalled as 
a genuine crisis of the CDM scenario fe.g. lElores fc Primackl 
1994; Moore 1991). Recent analyses show that observational 
data may actually be consistent with steeper profiles once 
the effects of inclination, non-circular orbits and triaxiality 
of the dark matter haloes are accounted for iHavashi et alJ 



^ Note that, in this region, numerical experiments are severeb 
hampered by two-body relaxation (sec e.g. Dicmand ct al._2004 
and discreteness effects iBinnevll2004li) . 
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Table 1. Radial density profiles of the Milky Way dark halo 
considered in the present work, according to the parametrization 
given by expression ISt . 



MW 


a 


/3 


7 


ro [kpc] 


po [GcV cm-3] 


ISO 


2 


2 





4.0 


1.655 


BE 


1 


3 


0.3 


10.2 


1.459 


NFW 


1 


3 


1 


16.7 


0.347 


M99 


1.5 


3 


1.5 


29.5 


0.0536 



2004. 120051). but t he controver sy is still unresolved (e.g. 



iGentile et al]|2004l : Ide 6 10)^ 2005). 

I n the Milky Way, it has been argued iBinnev fc Evan^ 
I2OOII . hereafter BE) that the microlensing optical depth to- 
wards the galactic centre reported by the MACHO collab- 
oration (Alc ock et a l. 2000) would be incompatible with a 
7 > 0.3 dark matter halo. A recent revision of the MA- 
CHO results yields a lower optical depth, in somewhat bet- 
ter agre ement with the v alues found by other experiments 
(see e.g. ISumi et al.ll2005l for a recent review) . Although a 
low optical depth could be consistent with steeper profiles, 
it is still unclear whether it would be compatible with the 
asymptotic slopes characteristic of numerical haloes (see e.g. 
[Birmcv 2004a). 

On the other han d, the existence of a b lack hole at the 
centre of our galaxy JCondolo fc Silklll999D and the adia- 
batic contraction of t he dark matter compon ent due to the 
presence of baryons jBlumenthal et allll98(fl are expected 
to increase the central dark matter density, leading to a very 
steep profile in the innermost regions. 

Most theoretical predictions of the gamma ray emission 
due to dark matter, including those performed in a super- 
symmetric framework (i.e. for heavy DM candidates), are 
based on 'cuspy' density profiles, since these maximize the 
expected flux. Given the present uncertainties, though, we 
have followed a completely different approach, trying to con- 
strain the shape of the Milky Way dark matter halo as inde- 
pendently as possible from the precise nature of dark matter 
particles. 

We therefore have considered four different models of 
the density profile of our galaxy, each one featuring a dif- 
ferent asymptotic slope at the centre: in order of decreasing 
'cuspiness', M99, NFW, BE and a non-singular isothermal 
sphere (hereafter ISO). Their corresponding parameters are 
summarized in Tabled and the density and cumulative mass 
profiles are depicted in Figure |5| In addition, we also con- 
sider a family of models in which a and (3 are fixed to 1 and 
3, respectively, while 7 is varied in uniform steps A7 = 0.05. 

In all cases, the normalization of the models, po, is set by 
imposing a local dark matter density p{rQ) — 0.3 GeV cm"'', 
with vq = 8.5 kpc. The characteristic radius ro has been 
chosen so that the virial radius and mass are JJvir ~ 260 kpc 
and Mvir « 10^^ Mq. Note that the ISO model can only 
approximately satisfy this condition. 

The characteristic velocity of dark matter particles is 
also a necessary ingredient in our model of the Milky Way, 
as the LDM annihilation cross-section associated to Z' ex- 
change explicitly depends on this quantity. In many studies, 
the rough estimate Vdm ~ 10~^c is assumed to be accu- 
rate enough, but this is definitely not true for our present 
analysis. The positron emission arising from the velocity- 
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Figure 2. Top panel: Densty profiles considered in the present 
work. All of them are described by equation 0, with values of the 
parameters given in Table Bottom panel: The corresponding 
cumulative mass profiles. 



dependent term of the cross-section will be sensitive to the 
product p^{r)a^{r). Since we are trying to constrain the 
shape of the density profile, it is extremely important to 
properly account for the radial variation of the velocity dis- 
persion. 

Indeed, the velocity dispersion profile of dark matter 
haloes is known to change systema tically with radius. As 
shown by iTavlor fc Navarrol (|20^), the phase-space den- 
sity profile found in N-body simulations follows an approx- 
imate power law over several ord ers of magn i tude in r. 
This result has been confirmed bv iRasia et alj (|20o3) and 
[Ascasibar ct alj i2004ft for haloes of very different mass. Fol- 
lowing the latter, the density and velocity dispersion profiles 
would be subject to the phenomenological relationship 



p{r) 



(r) 



10 



1.46±0.04 



V3 

vir 



-(—) 

I UviJ 



(4) 



where i?vir is the virial radius of the halo and we de- 
fine l/vir = GMvir/RviT- The cosmology-related quantities 
Pc = 1.36 X 10" Mq Mpc"^ and Mvir ~ ^^lOOpc-Rvir are 
the critical density and virial mass corresponding to a 'con- 
cordance' ACDM universe {Qm = 0.3, = 0.7, h = 0.7). 
Substituting i?vir ~ 260 kpc for the Milky Way, we obtain 
Mvir ~ 10^^ Mq and Kir « 130 km s"^ The velocity dis- 
persion profile of our galaxy is thus expected to vary with 
radius as 



6 X 10" 



p{r) 



1 GeV cm~^ V 1 kpc 



-1 2/3 



(5) 



This formula has been used by iBoehm fc AscasibaJ 
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r [kpc] 

Figure 3. Velocity dispersion profile obtained from the Jeans 
equation for each of the models in Tabled Grey line on the top 
panel shows expression O . 



lj2004) to estimate the characteristic velocity of dark matter 
particles. The radial dependence of the velocity dispersion 
significantly reduces the emission from the boson-exchange 
channel, being roughly equivalent to an effective flattening 
of the density profile. 

However, expression (|KJ is a mere fit to the coarse- 
grained phase-space density profiles found in numerical sim- 
ulations, and thus its validity has only been tested for 'cuspy' 
density profiles. A more self-consistent approach is followed 
in the present study, where we derive the velocity dispersion 
profiles from the spherically-symmetric Jeans equation. 



GM(r) 



p{r) dr ' '"^ ^ J. j,2 ' 

where Ura is the radial component of the velocity and 
the anisotropy parameter /3(r) is defined as /3(r) = 1 — 
{^e(''))/('''ra (''))■ Assuming no radial infall, an isotropic ve- 
locity ellipsoid and vanishing velocity dispersion at infinity. 



a (r) 



3 
p{r) 



p{r[ 



GM{r) 



(7) 



The corresponding profiles are plotted in Figure |21 As 
can be readily seen, equation ||KJ provides a fair approxi- 
mation for NFW and M99 models, but it is certainly not 
adequate for shallower density profiles. 



4 COMPARISON WITH SPI DATA 

According to the model outlined so far, the rate at which 
new positrons are created is 



where F corresponds to the decay rate Fd for decaying dark 
matter and to the annihilation rate Fa = (a + hv\^)nd-m* 
for annihilating LDM particles, with ndm* being the number 
density of dark matter anti-particles. 

These positrons will be relativistic at the moment of 
their creation (-5^+ ~ ni^^c?). However, they can efficiently 
lose their energy through coUisional ionization or excitation 
in neutral Hydrogen and by interaction with plasma waves 
in ionized interstellar medium. We will make the approx- 
imation that they can only travel a short distance before 
becoming non-relativistic and annihilate with an electron at 
rest. Such an assumption is not very realistic for the outer 
regions of the galaxy, but it is perfectly reasonable for the 
galactic bulge, where most of the observed emission comes 

from. 

Both OSSE ("Kinzer et al."200l'), TGRS jHarris et alJ 

Il998h and SPI ( Churazov ct al. 2005) measurements indi- 
cate that approximately 93 per cent of the positrons annihi- 
late through positronium formation. In this channel, 3/4 of 
the annihilations take place in the orthopositronium state, 
yielding 3 photons with E < 511 keV each, while the remain- 
ing 1/4 annihilate in the parapositronium state, producing 
2 photons with E = 511 keV. The remaining 7 per cent that 
do not form a positronium annihilate directly into 2 photons 
with E = 511 keV. Consequently, the total number of 511 
keV photons produced per unit time would be given by 



n-, = 2 (0.07 + 0.93/4) n„+ = 0.605 n„+ . 



(9) 



The predicted intensity distribution for any particular 
model of the Milky Way dark halo can be thus computed as 
the integral along the line of sight, as a function of galactic 
longitude I and latitude b, of the emissivity n-,{r), 



n-y{r) ds, 



(10) 



where the spatial dependence arises through the radial den- 
sity and relative velocity profiles Udmir) = ndm* (r) = 



/5dm(r-)/(2mdm) and vl^ 
at the earth is simply 



I{l,b) dfi. 



(r) « a (r). The total photon flux 



(11) 



(8) 



Although the resulting sky map will obviously preserve 
the spherical symmetry of our models, the morphology of 
the emission (more specifically, its concentration) depends 
on the shape of the Milky Way halo as well as on the DM 
annihilation cross-section or decay rate^. Observations are 
mostly sensitive to the details of the central part, where the 
intensity of the 511 keV line is highest. Dark matter in the 
outermost regions can make a significant contribution to the 
total fiux, but the intensity of the emission is so low that it is 
difficult to discriminate from the instrumental background. 

In order to compare with observational data, intensity 
maps I{1, b) have been computed for |/| < 60° and jbj < 50°. 
The total flux within this area is denoted by <[>tot . However, 
a fairer comparison with the flux measured by the satellite is 
given by the central 33° (~ 1 steradian) . We shall quote this 



■* Note that the precise distribution of baryons and their chem- 
ical composition should also be included in a realistic model of 



positron propagation. 
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Table 2. Theoretical photon fluxes (in cm ^ s ^) expected for 
different halo profiles and DM type. 



DM 


MW 






Td 


ISO 


0.0459 


0.0201 r26 




BE 


0.0439 


0.0195 r26 




NFW 


0.0478 


0.0232 r26 




M99 


0.0512 r26 


0.0262 r26 


a 


ISO 


5.13 026 


3.175 026 




BE 


5.21 a26 


3.325 a26 




NFW 


9.52 a26 


7.487 a26 




M99 


26.6 026 


24.5 a26 


b 


ISO 


1.92 X lO"'' fe26 


1.12 X 10"^ 626 




BE 


1.87 X 10-s 626 


1.11 X 10-'' 626 




NFW 


2.93 X lO"'' 626 


2.10 X 10-*^ 626 




M99 


5.90 X 10-^ 626 


4.97 X 10-6 626 



Table 3. Results of the model-fltting analysis. Fluxes in units 
of 10-3 cm-2 s-i. 



DM 


MW 


*tot 


*con 


MLR 




ISO 


6.82 ± 0.58 


2.95 ±0.25 


135.2 




BE 


7.23 ± 0.57 


3.18 ±0.25 


167.3 




NFW 


7.36 ± 0.46 


3.53 ±0.22 


261.2 




M99 


6.86 ±0.37 


3.48 ±0.19 


332.0 


a 


ISO 


5.55 ± 0.33 


3.40 ± 0.20 


282.8 




BE 


4.98 ± 0.27 


3.16 ±0.17 


353.6 




NFW 


2.49 ±0.11 


1.95 ±0.09 


459.9 




M99 


0.83 ± 0.04 


0.76 ± 0.04 


339.2 


b 


ISO 


6.00 ±0.38 


3.46 ±0.22 


258.3 




BE 


5.76 ± 0.32 


3.40 ±0.19 


305.7 




NFW 


3.61 ± 0.18 


2.57 ±0.13 


422.4 




M99 


1.57 ±0.07 


1.32 ±0.06 


430.0 



flux as $ccn. The values of $tot and <l?con expected for each 
combination of dark matter type and radial density profile 
are given in Table |5| 

Our analysis has been performed on the December 10, 
2004 public INTEGRAL data release, which consists of ^309 
days of observations. In order to reduce systematic uncer- 
tainties in the analysis, we exclude observation periods with 
strong instrumental background fluctuations^. The total ef- 
fective exposure time after cleaning is 15.3 Ms. The expo- 
sure is quite uniform in the central regions of our Galaxy 
< 50° and |6| < 30°). 

We use a maximum likelihood algorithm to compare the 
theoretical sky maps with the INTEGRAL/SPI data. This 
method has already been applied to SPI data to characterize 
the morphology of the annihilation. A detailed description 
can be found in lKnodlseder et al.lll2005D . 

Briefly, the normalization of each theoretical model 
is fitted to reproduce the measured rate in the 507.5 — 
514.5 keV energy range, taking into account an instrumen- 
tal background model, the pointings history and the spa- 
tial and energy response functions of SPI. Normalized maps 
have been convolved with the response function, providing 
the expected number of counts in each detector as a func- 
tion of the pointing periods. We then find the intensity that 
maximizes the log likelihood. We subtract from this log like- 
lihood Li the log likelihood Lo that is calculated under the 
hypothesis that there is no 511 keV source. Multiplication 
by a factor of 2 provides the maximum log-likelihood ratio, 
MLR = 2(Li-Lo). 

Results of the model-fitting procedure are presented in 
Table |21 and Table 3] As in the theoretical models, <E>tot is 
the total flux of the map, integrated over the whole solid 
angle, while "l>ccn is restricted to an aperture of 33°. When 
comparing two models, the one with the largest MLR can 
be said to explain the data better than the other, although 
differences Amlr < 10 are, in principle, not very significant. 

Strictly speaking, the posterior probability distribution 
for the flux, given the data and a particular model, is given 
by: 



^ These background variations are generally due to solar flares 
or exit and entry of the observatory in radiation belts. 



Table 4. MLR and total flux (normalized to lO-'^ cm-^ s-^) 
for the combined models ifiala + (pblbj where la and are the 
intensity distributions associated with a and 6-terms, respectively. 



DM 


<j)a 


06 


*con 


MLR 


ISO 


28.07 ±2.68 


-26.35 ±2.85 


1.72 


368.1 


BE 


15.11 ± 1.32 


-13.95 ± 1.53 


1.16 


437.1 


NFW 


2.42 ± 0.39 


-0.67 ±0.53 


1.75 


461.4 


M99 


-1.00 ±0.16 


2.82 ±0.25 


1.82 


464.8 



P{^\D,M) d$ = 



P{D\M,<i>)P{<i>\M) d$ 



(12) 



P(D\M) 

where D is the data, M is the model and <E> is the value of the 
flux, which can be regarded as a vector in the two-parameter 
models. P{D\M, <I>) is the likelihood of the data given M 
and P($|M) d<I> is the prior probability distribution of 
the flux, and P{D\M) is the Bayesian evidence of model M, 



P{D\M) 



P{D\M, ^)P(^\M) d<I>. 



(13) 



The Bayesian evidence provides a measure of how well 
the model explains the data. More precisely, model Mi 
can be considered P{D\M\) / P{D\Mo) times more likely 
than model Mq. Since the likehhood P{D\M, <I>) is strongly 
peaked around the best-fitting value of $, the evidence is 
not very sensitive to the precise shape assumed for the 
prior. Morevoer, for practical purposes we can approximate 
P($|_D,M) by a Gaussian centered around the best-fitting 
$ with a ~ A$, so that the evidence can be computed as 

P{D\M) oc e^^'' V2^A$ (14) 

for the models with one parameter, and 

P{D\M) oc e^^"- 27rA(/>aA0i, (15) 

for the composite models in Table 2] 

It is evident from the figures in Table El that we can 
neglect the effect of A$ when comparing two models. In 
other words, a fair approximation to the Bayes' factor can 
be obtained by simply comparing the likelihood ratios, 

P{D\Mi) 



log ■ 



(MLR)i - (MLR)o. 



(16) 



P{D\Mo) 

The model that best fits the INTEGRAL/SPI is a 
NFW profile with constant cross-section. This model closely 
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ressembles the halo model H' of iKnodlseder et al.l i2005l) . 

with similar values of the best-fitting flux and MLR. All 
other models in Table have Amlr > 20 with respect to 
the g-term NFW or an y of the bulge and halo models in 
IKnodlseder et all ll2005l) that satisfactorily explain the ob- 
served morphology of the 511 keV emission. In frequentist 
terms, all the other density profiles and dark matter types 
listed in Table|3|can be rejected at the (1 — e"^") confidence 



Substituting 1191 1 into 1171 1 yields 



level, or ~ \/20a. 

A rigorous Bayesian comparison between the one- and 
two-parameter models would require an extensive discussion 
of the priors. However, we note that linear combinations of 
the a-term and &-term models (c.f. Table^J lead to marginal 
improvements of the MLR, but they provide non-physical 
solutions where one of the fiuxes is negative. When the fit is 
constrained to non-negative values, zero fluxes are obtained 
for the negative coefficients, and the results quoted in Ta- 
ble |3 are reproduced for the other component. We therefore 
conclude that two-parameter models fail to provide a better 
explanation of the observational data. 



5 RESULTS 

5.1 The nature of dark matter 

Comparing the best-fitting flux in Tablej^for a NFW profile 
with an a-term with our theoretical prediction in Table |21 
we obtain 



026 = (2.6 ±0.12) X 10 



-4 2 

^McV- 



(17) 



This result is valid for any kind of candidates, as long 
as the annihilation cross-section is parameterized as we have 
done in the present work and the galactic positrons are pro- 
duced as primaries. 

Equation 1171 becomes comparable with Q when 
mdm — 100 MeV. This confirms our previous conclusions 
from the gamma ray constraints, i.e. that the LDM scenario 
requires both an a and &-term when rridm < 100 MeV. 

To explain the 511 keV emission line, candidates heav- 
ier than 100 MeV would require a (velocity-independent) 
annihilation cross-section into a pair electron-positron that 
is well above the relic density requirement. Heavy candidates 
must thus produce positrons as secondaries to be a possible 
solution. However, one would expect an overproduction of 
gamma rays in the Milky Way in such case, unless there ex- 
isted a channel that lead to a large production of positrons 
and a low production of gamma rays. This is why we advo- 
cate for LDM particles with mdm < 100 MeV in order to 
explain the origin of galactic positrons. 

For a scalar particle, the F-exchange channel yields 



c 

47V 



WdnJ WdiJ 



I 2 



„2 2 



(18) 



where mf ^ mdm is the mass of the fermions in the final 
state, mF is the mass of the heavy fermion that is exchanged 
during the annihilation and the quantities c; and Cr corre- 
spond to the couplings. 

Neglecting m/ compared to mdm, we flnd 



4 TV fn?p 



a26 



9.32x10^ cfcl ( 



rap 



100 GeV 



(19) 



X 10'' 



C; Cr 



100 GeV mMov' 
where the two couplings ci and Cr are expected to be lower 
than unity (a few units at most). 

For c; « Cr ~ 1, one obtains mp ^ 6 — 0(600) TeV 
(for 1 < mMcV ^ 100). This is obviously out of reach 
for past and forthcoming colliders. However, for smaller 
but more realistic values of the couplings, rap could be 
within the range of next colliders. According to equation 
l|2nji a mass mp £ [lOOGeV, O(TeV)] would correspond to 
ci Cr € 1.67 X [10"'*, 10""^] mMcV. Those are not particularly 
small couplings (especially for mMcV ~ 100 MeV) so there 
might be a signature of the LDM scenario in the next gen- 
eration of accelerators, and in particular at LHC if the F 
particles associated with quarks {Fq) are not too heavy and 
their couplings not too small. Note that smaller couplings 
e — Fe — dm (leading to mFe < 100 GeV) are forbidden due 
to the absence of signal in past colliders (notably LEP). 

So far we have considered scalar LDM particles. How- 
ever, fermionic candi dates (either Dirac or Majorana) have 
also been suggested jBoehm fc Favel]l2004ll . The difference 
with respect to the scalar case would be that, for a Dirac 
fermion, the cross-section would be given by 



aD ' 



■77^ (Ci + Cr) 

62 TV 



instead of H19|l . while for Majorana particles, 



— (Ci Cr) 

Stt m 



(21) 



(22) 



where F now denotes a scalar. 

Using equation l)17|l , we obtain that the mass mp should 

be 



mp 



0.145 y/{c 



100 GeV 

for Dirac dark matter particles and 

0.206 Cl Cr 



(23) 



(24) 



100 GeV 
for Majorana candidates. 

Assuming realistic values for the couplings (i.e. ci^r < 
0(1)), we find that the F mass is much smaller than 100 GeV 
regardless of whether LDM is a Dirac or Majorana fermion. 
Since the presence of charged particles much lighter than 
~ 100 GeV has been excluded by LEP data, one readily 
sees that fermionic LDM particles cannot explain the 511 
keV line emission unless one considers couplings at the edge 
of perturbativity. 

Our results also indicate that a Z' cross-section can- 
not explain the observed 511 keV emission on its own. A 
cross-section strictly proportional to can be ruled out by 
Amlr ^ 29.9. The best-fitting fiuxes obtained for a linear 
combination of the models with an a- and a 6-term are given 
in Table 2] where two free parameters have been adjusted 
to match the SPI data. Although the MLR is somewhat im- 
proved, a non-physical solution (with one of the factors be- 
ing negative) is always obtained, suggesting that the boson- 
exchange channel plays only a minor role within the Milky 
Way halo. 

The existence of scalar DM coupled to heavy (fermionic) 
particles is thus required for our model in order to explain 
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Figure 4. Amlr = MLR(7) - MLR(l) for different central 
asymptotic slopes 7. Continuous line plots the best-fitting ninth- 
order polynomial. 



the 511 keV l ine, a t variance with the results reported by 
iBoehm et alJ (|2004 ) where the F were th ought to be fac - 
ultative, but in agreement with lBoehm fc A scasibar! i2004) . 
The reason for the discrepancy resides in the more accurate 
description of the velocity dispersion profile of the Milky 
Way. 

For the density profiles considered in the present work, 
decaying dark matter is completely incompatible with the 
observed morphology of the 511 keV emission. Even for the 
best-fitting model, M99, decaying dark matter of any sort 
can be ruled out by Amlr = 127.9. It would be possible, 
though, that a steeper dark matter halo (e.g. as predicted 
from adiabatic contractio n) may provide a better match to 
the observations (see e.g. lPrada et al.ll2004ll . 



5.2 The Milky Way density profile 

All of the configurations in Table|3|display Amlr > 20 when 
compared to our best-fitting model, namely a NFW halo 
of annihilating dark matter with approximately constant 
cross-section. We will now assume that this is indeed the 
nature of dark matter particles. If most galactic positrons 
did arise from DM annihilations, the spatial distribution of 
the 511 keV emission line would provide an extremely sen- 
sitive probe of the shape of Milky Way dark halo. 

We plot in Figure 2] the MLR of a series of models in 
which the parameters a and j3 in expression have been 
set to Q = 1 and /3 = 3, whereas the inner slope 7 varies from 
7 = 0.1 to 7 = 1.6 in uniform steps A7 = 0.05. We have 
normalized the MLR to the particular case 7 = 1, which cor- 
responds to the NFW profile. Fitting the data points with a 
ninth-order polynomial, the optimal value of the logarithmic 
slope is found to be 7 = 1.03 ± 0.04, where the errors have 
been estimated by equating Amlr = 1. 

This result is perfectly compatible with 7 = 1, but any 
of the other profiles suggested in the literature would be 
extremely hard to reconcile with the INTEGRAL /SPI data. 
This is again at odds with iBoehm et al.l i2004) . where a 
shallower profile with 7 ~ 0.6 was favoured, based on a 
coarser comparison between the theoretical predictions and 
the observed fiux and extension of the emission. 



Finally, we would like to stress that the constraint we 
obtain for the inner asymptotic slope of the density profile is 
so tight that, if the Milky Way dark halo was found to follow 
a different shape by some independent means, the possibil- 
ity that dark matter annihilations were the main source of 
galactic positrons would seem rather unlikely. Systematic ef- 
fects (see Section f6. in would in general tend to yield values 
of 7 below the real one, so our estimate should be regarded, 
to a certain extent, as a lower limit. If DM is responsible 
for the 511 keV emission, 7 > 1. If 7 < 1, galactic positrons 
must come from a different physical process. 



6 DISCUSSION 

6.1 Systematic effects 

The morphology of the galactic 511 keV line emission pro- 
vides a wealth of information on both the nature of dark 
matter and the shape of the Milky Way dark halo. We have 
shown in the previous section that several constraints can 
be derived for the parameters that characterize dark mat- 
ter particles (with special emphasis on the LDM scenario), 
as well as the inner logarithmic slope of the radial density 
profile. 

One should keep in mind that our analysis is based on 
several simplifying assumptions. Relaxing one (or each) of 
them would have a different systematic effect on our results: 

First, we have neglected any astrophysical contribution 
to the galactic positron budget. Emission from sources other 
than dark matter would lower our estimate of a26»7i^^y, 
equation 11711 . by a factor proportional to the fraction of 
DM-related positrons. The effect on the density profile de- 
rived for the Milky Way dark halo depends on the spatial 
distribution of the other sources. Our estimate of the inner 
slope 7 would be biased low if the latter was flat, and high 
if the sources were concentrated near the centre. 

Second, the fact that positrons may travel a certain dis- 
tance before losing their energy and annihilating would flat- 
ten the expected emission. Therefore, a steeper density pro- 
file would be required in order to fit the observations (i.e. 
our estimate of 7 would be biased low). The best-fitting flux 
would be somewhat lower, so the real value of 026 ''tt-iacv 
would again be lower than (1171 . 

Third, our model of the Milky Way is overly simplis- 
tic in several respects. On one hand, it is well known that 
the dark matter haloes found in numerical simulations dis- 
play a significant degree of triaxiality (Jing &: Suto 2002j), 
although the incl usion of gas cooling ten ds to yield more 
spherical haloes jKazantzidis et all 12004) . In our galaxy, 
observations of the Sagittarius tidal stre am have been in- 
terpretcd as favouring a spherical halo ilbata et al."200l|; 
Majcws ki ct al. 2Qo3), but recent analy ses also suggest both 
oblat e iMarti'nez-Delgado et alJl2004l) and prolate (IHelmil 
I2OO4I) shapes. Although the INTEGRAL/SPI data is con- 
sistent with spherical symmetry, higher-resolution observa- 
tions would be needed in order to quantitatively address this 
issue. 

In addition, some substructure is expected to be present 
in the dark matter halo of our galaxy, both in real and 
phase space (i.e. the six-dimensional space of positions and 
velocities). In real space, dark matter clumps would tend 
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to boost th e expected emission d ue to the increase in lo- 
cal density jBere^strom et alJll999l) . It has been argued (e.g. 
iHooner et alJl2004h that the 511 keV line could actually be 
detected not only from the galactic centre, but also from the 
nearest dwarf spheroidals. On the other hand, structures in 
phase space (such as tidal streams) would not have any effect 
on the emission through the a-channel, but they may have 
a significant impact both on direct detection experiments 
iHelmi et alJ |2002) and on the emission arising from the 
6-term of the cross-section, mostly because of the increase 
in dark matter velocity with respect to the local velocity 
dispersion. Although it does not seem likely, for 626 ~ 10, 
that structures in velocity space yield a detectable signa- 
ture, there might be a mild enhancement in the signal from 
local dwarfs, particularly near their pericentre. 

A more realistic model of the density and velocity dis- 
tribution of dark matter particles within the Milky Way 
halo would be given by the results of N-body simula- 
tions. Althou gh extremely promising, this approach (see e.g. 
IStoehr et al] l2003) must face the problem of numerical con- 
vergence (i.e. lack of resolution) in order to provide conclu- 
sive results. 



with the gamma-ray data. To make this point, they wrote 
the FSR cross-section as 



dob 
dE. 



, — CI" ann X -Rcorr (-£/) 



and used R^. 



In 



fe) 



1 



(25) 



with 



s = 4mdm and s' = 4m dm (mdm — E-y), as computed 
by iBerends fc Bohml il987ft for the radiative correction to 
e"'"e~ — > fj,'^ iJ.~ . 

To our knowledge, it has never been proven that any 
radiative cross-section can be written in the same way as 
e'^e^ — > in the hard photon limit, so it might be pos- 

sible that equation l|25|l does not provide a completely accu- 
rate description of the radiative corrections to the annihila- 
tion process. Nonetheless, a detailed computation would be 
quite delicate, and it certainly lies well beyond the scope of 
the present work. We will therefore estimate t he gamma ray 
flux according to the procedure followed bv lBeacom et all 
(2005). Assuming the same Rcott{E), the intensity of inter- 
nal bremsstrahlung gamma rays at the earth is related to 
the 511 keV sky maps by 



6.2 Compatibiliy with low-energy gamma rays 

DM annihilations produce charged particles which, in turn, 
produce low-energy gammas. But low-energy photons are 
also expected through the final state radiation (FSR) mech- 
anism or DM annihilations into two monoenergetic photons 
(with an energy E — mdmC^). In principle, the observa- 
tion of such line would be an unambiguous signature of the 
LDM scenario. However, the corresponding cross-section re- 
lies on a box diagram that involves four powers of the electric 
charge e, and therefore it is expected to be suppressed com- 
pared to the FSR mechanism in which either the electron, 
the positron or the charged particle that is exchanged during 
the annihilation emits a photon. 

A rough estim ate of the gamma ra y flux generated by 
FSR was given in iBoehm et alJ i2004ft . In an attempt to 
build a model that would surely satisfy astrophysical con- 
straints, the authors assumed that the FSR cross-section 
was as large as the annihilation cross-section. The electric 
charge (which is expected to reduce the FSR cross-section by 
a factor e^) was therefore deliberately omitted in order to ob- 
tain the maximal flux of low-energy gamma rays that LDM 
could produc e . Com paring this estimate with observations, 
iBoehm et alJ (|200^ deduced that particles lighter than 100 
MeV were excluded unless their annihilation cross-section 
was suppressed in our galaxy with respect to its value in the 
primordial universe, which led the authors to postulate the 
existence of a velocity-dependent annihilation cross-section. 

The required suppression factor turns out to be consis- 
tent with that obtained from the analysis of the 511 keV 
line. More precisely, a cross-section as large as equation I17II 
does not overproduce gamma rays if the DM particles are 
lighter than 30 MeV. It does, in principle, if DM is heavier, 
but given the drastic assumption on the FSR cross-section, 
it was considered that even 100 MeV particles would both 
explain the 511 keV line and be compatible with gamma-ray 
observations. 

In a recent paper, iBeacom et al.l (|200^ claimed that 
particles heavier than 20 MeV were actually not consistent 



d/b 



-^511 



1 + (s'/sf 



In 



(26) 



Comparison with COMPTEL and EGRET data only 
requires the computation of the average intensity over the 
appropriate region of the sky. Using the model that best flts 
the 511 keV line, i.e. a NFW profile with constant cross- 
section given by equation 1171 . we obtain 



|<30 ,|b|<5 



^ 0.005 cm~^ s~^ sr"\ (27) 



Even if we restrict to the innermost 5° in both longitude 
and latitude, our expectation for the average intensity 



: 0.018 cm" 



sr 



(28) 



;|<5Mbl<5 

is som ewhat lower than the estimate used bv lBeacom et all 

The solid angle considered and the precise shape of the 
energy dependence of the observational data have important 
consequences on the inferred constraints. Figure|^shows the 
limits we obtain for the mass of the dark matter candidate. 
The factor of ~ 4 difference between equations 1271 and 1281 1 
translates into a factor of 3 for the maximum dark matter 
mass, i.e. rridm < [30 — 100] MeV. This con straint is simi- 
lar (though less severe) to that obtained bv lBeacom et all 



J2005I). and it al s o agr ees roughly with the results reported 
bv lBoehm et alJ (|22q3) despite the fact that much less than 



2 photons are predicted by integrating equation 1261 over 
the photon energy (e.g. between 1 and 100 MeV). However, 
the constrain provided by the 'spectrum' of the observed 
gamma-ray radiation is much more stringent than its total 
flux (or, more prec isely, the flux above ^ mdmC^, see 

iBoehm et 'ai]l2004^ . 

In any case, we would also like to emphasize that a 
rigorous comparison with COMPTEL and EGRET data, 
analogous to the analysis of the 511 keV line presented here, 
would be necessary in order to provide accurate constrains 
on the mass range allowed for the dark matter particles. 
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E [MeV] 

Figure 5. Comparis on w ith COMPTEL (circles) and EGRET 
(diamonds) data (.Stro ng et al -':200o|). Lines plot the gamma ray 
flux expected for ra^^^ = 30 and 100 MeV, according to expres- 
sion l26l Solid and dashed lines correspond to the regions j6| < 5°, 
\l\ < 30° 03 and |b| < 5°, \l\ < 5° FM . 



6.3 Tests of the LDM model in particle physics 
experiments 

We mentioned in a previous section that the F particles 
could be detectable in the next colliders. The F particles 
could be detected through their two-body decay Fp~ —> p~ + 
where denotes the missing energy associated with Dark 
Matter and p~ the Standard Model particle associated with 



Fp (for example Fe 



e + ^) • The couplings are large 



enough to allow for the decay within the detector, although 
the production cross-section of _F particles {p'^p~ F~ F^) 
may be too small to yield a visible signature. This cross- 
section can reach, however, a few pb for m^p ~ 100 GeV, 
T7^dm ~ 100 MeV and couplings Fp — p — dm of the order of 
12011 . The presence of Fq particles could then be detected at 
LHC through the modification of the total hadronic cross- 
section (if either of the couplings ci or Cr remains relatively 
large) . 

The Fe decay should be quite similar to the (supersym- 
metric) decay of a chargino into an electron and a sneutrino. 
However, in our case the mass of the DM particle is much 
smaller than the sneutrino mass. Hence, we expect a fair 
repartition of energy between the electron and the missing 
energy in the LDM scenario. This should be quite different 
in a supersymmetric framework, albeit more precise esti- 
mates would depend on the difference of mass between the 
sneutrino and the chargino (and therefore on the model con- 
sidered). In any case, a possible way to discriminate between 
our scenario and supersymmetric particles could be the ab- 
sence of other signatures. For example, unlike the chargino 
which is expected to also decay into a selectron and a neu- 
trino, the Fe particles are assumed to leave only one main 
signature (i.e. the decay into an electron and missing en- 
ergy). The precise relationship between the Fe mass and the 



couplings 12011 should also help in discriminating between 
the LDM and supersymmetric scenarios. 

Dark matter production through p'^p~ — > dm dm could 
in principle be detectable in the initial state radiation pro- 
cess (where a photon is emitted by the particles of the initial 
state). However, the associated cross-section is expected to 
be smaller than a few fb even for mdm ~ 100 MeV, and is 
therefore less interesting. 

There are two other tests that could be really crucial for 
our purpose. One is based on the so-called NuTeV anomaly 
and the other one on the value of the fine structure constant 

Q. 

NuTeV is an experiment which measured the ratio 



R 



neutral currents 2 
charged currents 



with gl = [{gr,,f + {gt.f]/4 and 5/;, = 2{n{fi.r) - 
Q{f) sind^) the left and right couplings of the Z boson to 
fermions. Through the measurement of this ratio, one can 
infer the value of the mixing angle sin 6^ and compare it to 
Standard Model expectations (obtained notably from LEP 
electroweak precision observables) . Surprisingly, NuTeV re- 
sult for sin^^ turns o ut to be slightly discrepant with the 
Standard Model value fMcFarlan d et alJl2003l: iNaples et alJ 
2003; McFarland & Moch 2003). 

There are several possible explanations to the NuTeV 
anomaly. In particular, isospin violation and/or strange sea 
cisymmetry as well as other effects such as electroweak cor- 
rections may reduce significantly the discrepancy. However, 
at present, the situation is still uncertain and there is an 
open window for new physics. The best exp lanation to that 
respect turns out to be a light gauge boson ijDavidson et alJ 
2002; Bochm 2004), similar to the one we introduce in the 
LDM scenario. If the NuTeV anomaly disappears or is sig- 
nificantly reduced, then one could set a limit on the Z' cou- 
plings to ordinary matter. On the other hand, more evidence 
in favour of a light Z' would certainly give more credit to 
the LDM scenario. 

The other t est concerns the existence of heavy particles 
ijBoehm fc Ascasibar 20o3l . We expect them to contribute 
to the electron anomalous magnetic moment as: (5ap,e ~ 



167r^ mp 



> 0. Using expression 1201 . we obtain 



5.41 X 10 



^12 '^"''dn 



MeV 



(29) 



It turns out that there is a small discrepancy between 
the theoretical value of (hereafter denotes ath) and its 
measurement (oe^p): 



ath) ~ (3.44 3.49) 10' 



(30) 



where the first immber is obtained from the positron g- 
2, while the second one is from the electron. We estimate 
1301 1 by using ath = /(a) with a — qqh, the fine structure 
constant as measured by Quantum Hall efi'ect (QH) experi- 
ments. There are other experiments aiming at measuring a, 
but QH experiment s seem the most precise at present (see 
e.g. lKinoshitalll996D . 

Usually, one assumes instead the validity of QED and 
imposes that ath = agED matches ae^p- One then gets a 
'theoretical' estimate of a that is extremely precise and in 
fact the most significant input in the a value given in the 
CODATA. The latter differs too much from agn for the dif- 
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ference to be explained by common extensions of the Stan- 
dard Model such as supersymmetry, but experiments mea- 
suring Q did not reach the sensitivity of gr — 2 experiments 
as yet. Therefore, there is the hope that the difference may 
eventually go away and that there is no new contribution 
other than QED. 

However, the discrepancy between aqn and a could 
originate from new physics. In particular, the introduction 
of heavy particles coupled to light scalars adds a new con- 
tribution in ath which is greater than expected in e.g. su- 
persymmetry, due to the DM mass scale. 

Including this new contribution iSa^), one obtains a ~ 
O.QH (or ath — aexp) if mdm ~ 6.4 MeV , using the aver - 
age value of ath computed with aqn (see lKinoshital[T99d) . 
Taking into account theoretical and experimental uncertain- 
ties, mdm ~ 3 — 9 MeV. For smaller DM masses, we obtain 
5a^ ^ Afle, while for larger DM masses 5a^ ^ Aa^. As ex- 
plained above, there is no direct measurement of a that is 
as precise as the gr — 2 as yet. Therefore, it is hard to exclude 
values above 7 MeV. However, this certainly places a very 
strong contraint and motivates further experiments measur- 
ing the value of the fine structure constant directly (and 
independently of QED). If these experiments find a perfect 
agreement with the value recommended in the CODATA, 
then the LDM scenario will have difficulties in explaining the 
511 keV line emission. If they found a discrepancy (whether 
it is positive or negative) then LDM will remain a serious 
candidate because it would be the sign of new physics. In 
particular, if the value oqh is confirmed, then the LDM 
scenario may reconcile the results from both g — 2 and a 
experiments, despite the difference of sensitivity. 

Measurements of the fine structure interval of ffe-like 
ions based on laser spectroscopy might also provide a very 
precise determination of a if theory gets as accurat e as the 
exper imental determinations (see e.g. Pachucki & S apirsteiiil 
l2003l . and references therein). However, this method assumes 
that there is no additional contribution from new physics, so 
it may not be suitable for answering whether the F particles 
exist or not. 

Taking the same couplings and the same mass mp for 
Fe as for F^, we obtain a very large contribution to the 
muon g — 2. Our prediction, in fact, exceeds the exper- 
imental value by a factor 2-3, which is itself larg er than 
the Standard Model prediction iBennett et aljl2004) . It was 
found Aa^ = {a^xp — ath) £ [1.6,2.7] 10~^. So, by using 
»71mgV ~ 6 — 7 and mpii = Smpe (or e.g. mpp. = 2mFe and 
smaller couplings to the muons), our prediction for the muon 
g — 2 becomes compatible with the experimental value. In 
fact, the LDM scenario would even explain the well-known 
discrepancy. Note that such a hierarchy exists in the Stan- 
dard Model and it is very realistic to assume that it exists 
also in any other extensions. 

Hence, the LDM scenario could in fact explain both the 
experimental values of the fine structure constant and the 
muon g — 2 for mdm 6 — 7 MeV. 



7 CONCLUSIONS 

In this paper, we use the intensity and morphology of the 
observed 511 keV line to put independent constraints on 
the nature of LDM and the shape of the Milky Way dark 



halo. Our main assumptions are that DM decays or anni- 
hilations are the main source of galactic positrons and that 
these positrons do not travel long distances before annihi- 
lating. Theoretical expectations for the flux distribution are 
computed for different DM models and galactic density pro- 
files. The sky maps thus obtained are then convolved with 
the SPI response function and used as a source for the IN- 
TEGRAL model-fitting analysis. 

We can rule out from a likelihood analysis the possibil- 
ity that decaying dark matter is responsible for the observed 
emission, unless the density profile of the Milky Way dark 
halo turns out to be extremely cuspy (with inner asymptotic 
slope 7 > 1.5). 

We can exclude fermionic LDM particles, because it 
would require the introduction of charged scalars lighter 
than 100 GeV, which should have already been detected 
in past colliders. As a result, LDM is likely to be a scalar. 

For annihilating scalar LDM, it is shown that the ex- 
change of a heavy fermion {F^) is required in order to fit 
the morphology of the 511 keV line, while the existence of 
a Z' boson would be necessary to satisfy the relic density 
criterion. Assuming a full spectrum and, most precisely, the 
existence of F particles associated with quarks {Fq), we no- 
tice that there might be a signature at the Large Hadron 
Collider, notably through the measure of the total hadronic 
cross-section and the two body decay of these Fq particles. 

However, the most promising signature of F particles 
turns out to be their contribution to the electron g ~2. The 
new contribution would make the measurements of the fine 
structure by the Quantum Hall experiment and the electron 
anomalous magnetic moment compatible for mdm ~ 6 — 
7 MeV, meaning that the value of a quoted in the CODATA 
(and used for many estimates) may not be the correct one. 
The dark matter mass could however be larger than ~ 6 
MeV. The existence of clumps, the fact that the dark halo 
is probably not perfectly spherical and the contribution to 
the positron population from astrophysical sources are all 
expected to decrease our estimate in equation 1171 therefore 
allowing for larger dark matter masses. 

Assuming the existence of this spectrum (and Ff_i parti- 
cles), we also find a non- negligible contribution to the muon 
g — 2. Both and could then explain the discrepancy be- 
tween the Standard Model predictions and the experimental 
values of the muon g — 2 and the fine structure constant. Al- 
ternatively those could provide a way to constrain the LDM 
scenario. 

Concerning the shape of the Milky Way dark halo, our 
results clearly indicate that dark matter particles can only 
explain the observed 511 keV emission if our galaxy features 
a cuspy density profile. For any annihilating DM candidate 
with constant cross-section, the best-fitting inner asymp- 
totic slope is found to be 7 = 1.03 ± 0.04. 

To sum up, we would like to stress the fact that the 511 
keV emission line provides extremely stringent constraints 
on the light dark matter parameters. Independent confir- 
mations are needed to prove that dark matter contributes 
to most of the galactic positrons. Such confirmation might 
come either from the lack of astrophysical sources and/or 
from detection in particle physics experiments. Observations 
of the density profile of the Milky Way have the possibility 
to rule out a dark-matter related origin of galactic positrons 
if the density profile of our galaxy is found to be shallow 



12 Y. Ascasibar et al. 



at the centre. Alternatively, the discovery of LDM particles 
would have a tremendous impact on the determination of 
the dark halo profile of the Milky Way. 
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